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Abstract 


When  an  inventory  of  multiple  items  is  used  to  support  the  completion  of  an 
overall  job  or  mission,  commonly  used  inventory  service  level  measures 
which  are  defined  in  terms  of  item  availability  are  not  appropriate.  This 
paper  develops  a  general  model  for  problems  of  this  type  in  which  item 
demands  are  also  interdependent.  A  heuristic  is  derived  for  determining 
stock  levels  that  are  guaranteed  to  satisfy  the  given  inventory  system 
performance  objective,  but  at  a  slightly  higher  cost  than  would  be  achieved 
by  the  true  optimum.  ( ^ _ 


This  paper  describes  work  in  progress  and  should  not  be  cited  without  the 
authors'  consent.  Comments  are  welcome. 
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1.  Introduction 

Service  levels  of  multi- item  inventory  systems  are  often  defined  in  terms 
of  general  performance  measures,  such  as  the  number  of  "jobs  completed"  or 
the  number  of  "end  products"  produced,  where  each  job  or  end  product 
requires  a  combination  of  the  basic  inventory  items  for  completion.  For 
example,  repairing  a  machine  in  the  field  may  require  several  parts  and 
tools,  any  one  of  which  is  sufficient  to  cause  an  incomplete  job  if  it  is 
not  carried  in  the  repair  service  kit.  In  this  context,  the  performance  of 
the  repair  kit  is  defined  as  the  fraction  of  jobs  for  which  it  allows 
completion,  rather  than  the  fraction  of  part  demands  that  it  is  able  to 
satisfy.  In  the  context  of  field  repair  systems,  several  stocking  analysis 
methodologies  have  been  developed  for  this  problem.  [Smith,  Chambers  and 
Shlifer(1980) ,  Mamer  and  Smith(1982) , (1985) ,  Graves(1982) ,  Hausman(1982) , 
Schaefer (198 3) ,  March  and  Scudder(1984)  ].  These  approaches  recognize  the 
interdependencies  between  the  item  demands  in  multi -item  systems,  but 
generally  restrict  item  stock  levels  to  be  0  or  1.  This  restriction, 
which  is  imposed  to  obtain  computationally  feasible  solution  methods, 
limits  the  models'  application  to  items  which  are  restocked  after  each  job 
or  to  items  such  as  tools,  whose  inventories  are  not  depleted  after  use. 

In  recent  articles,  Baker,  Magazine  and  Nuttle(1986)  and  Baker(1985) 
analyze  a  structurally  similar  problem  in  the  manufacturing  context.  Their 
papers  consider  a  production  system  in  which  components  are  obtained  in 
advance  for  uncertain  demand.  Later  these  components  are  assembled  on 
order  into  final  products.  The  service  level  of  the  system  is  defltibd  as 
the  probability  of  meeting  all  final  product  orders,  while  inventory 
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Investment  is  determined  by  the  stock  levels  of  the  component  items.  These 
euthors  analyze  examples  to  illustrate  the  effect  of  component  commonality 
on  the  cost  of  meeting  specified  final  product  service  level  targets. 
Despite  the  structural  similarity  between  the  manufacturing  and  field 
service  problems,  the  previously  discussed  repair  kit  optimization 
techniques  do  not  apply  to  the  manufacturing  problem  as  stated,  because  the 
item  stock  levels  are  greater  than  one  and  the  component  inventories  are 
gradually  depleted  as  end  product  orders  are  filled.  Baker,  Magazine  and 
Nuttle  note  the  complexity  of  solving  general  problems  of  this  type  and  the 
need  for  more  research  in  multi-item,  multi-level  inventory  systems. 

In  this  paper,  we  propose  a  heuristic  for  extending  the  previously 
developed  job  completion  inventory  models  to  the  case  of  general  stock 
levels  with  inventory  depletion.  Because  of  the  complexity  of  the  problem, 
we  believe  that  a  direct  solution  even  for  moderately  large  numbers  of 
items  and  jobs  (or  final  products  in  the  case  of  manufacturing)  is  not 
computationally  feasible.  Instead  we  develop  an  objective  function  that 
serves  as  a  lover  bound  function  for  the  service  level  that  is  achieved  for 
a  given  stock  level.  The  minimum  cost  stock  levels  can  then  be  found  that 
will  bring  the  lover  bound  up  to  a  specified  service  level  target.  This 
approach  guarantees  that  the  service  level  objective  will  be  met  but, 
because  of  the  use  of  the  lover  bound  will  tend  to  result  in  additional 
inventory  cost.  Alternatively,  for  a  fixed  budget  constraint,  we  can 
maximize  the  service  level  bound.  This  second  case  is  somewhat  more 
powerful,  because  the  resulting  stock  level  choice  can  then  be  plugged  into 
the  actual  objective  function,  which  could  be  evaluated  by  Monte  Carlo. 

This  second  case  is  illustrated  by  an  example,  and  other  Intermediate 
results  are  illustrated  as  well. 

2.  Model  Specifications 

Key  data  for  the  problem  we  consider  is  contained  in  the  "job  matrix"  J, 
which  is  defined  as  follows: 

Jjj  -  the  number  of  items  of  type  i  required  to  complete  a  job 
(or  satisfy  an  order)  of  type  j. 

The  structure  of  a  job  matrix  is  illustrated  in  Figure  2.1  on  the  next 
page,  which  will  serve  an  example  for  subsequent  illustrative 
calculations.  In  general,  we  will  assume  there  are  n  item  types  (  i  - 


l,...,n)  and  a  job  types  (j  -  l,...,a). 


Jobs  or  end  product  orders  ere  sssuaed  to  fora  e  sequence  of  independent 
events,  whose  intererrival  times  have  aean  1/A.  That  is, 

A  -  rate  per  unit  tiae  at  which  requests  or  jobs  arrive. 
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Figure  2.1.  A  Job  Matrix  J 


Job  types  on  successive  arrivals  are  assuaed  to  be  independent.  The 
marginal  probabilities  for  job  types  are  defined  as  follows: 

pj  -  Plan  arriving  job  or  request  is  of  type  J). 

In  soae  cases,  we  will  assume  that  the  job  arrivals  are  also  Poisson 
processes . 

The  expected  deaand  rates  for  individual  items  are  determined  by: 

*i“AJJijPj*  (21> 

where 


Ai  »  average  deaand  rate  per  unit  time  for  item  i,  i-l,...n. 
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For  any  fixed  time  t,  let  us  define  the  random  variables: 

Nj(t)  -  number  of  requests  (jobs)  of  type  j  arriving  by  time  t, 

j  -  1 . b- 

X^(t)  -  number  of  items  of  type  i  requested  by  time  t,  i  -  l,...,n. 
Letting 

X(t)  -  Xx(t) . Xn(t) 

N(t)  -  N^t) . NB(t), 

we  clearly  have  the  matrix  relationship 

X(t)-JN(t).  (2.2) 

For  a  given  vector  s  -  s^ . sn  of  stock  levels, 

&£  -  Initial  stock  level  for  item  type  i, 

we  also  have 

P{all  Jobs  completed  up  to  time  t)  -  P(X(t)  S  s).  (2.3) 

Equation  (2.3)  thus  defines  one  performance  measure  for  the  multi -item 
system  with  stock  level  s;  namely  the  probability  that  a  time  t  will 
elapse  with  no  unsatisfied  requests. 

Another  performance  criterion  of  interest  is: 

E[ time  before  a  stockout  occurs]  -  r 


-  |  [1  -  P(X(t)  s  s) ]dt.  (2.4) 

The  evaluation  of  the  inventory  system's  performance  can  be  done  exactly  if 
the  probability  F(X(t)  5  s)  In  Equation  (2.3)  can  be  calculated. 

This  calculation  is  more  difficult  than  it  may  first  appear,  however.  For 
example,  suppose  that  the  job  arrivals  are  independent  Poisson  processes. 
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Then  N(t)  and  X(t)  ere  clearly  vector  valued  Markov  processes,  which 
have  renewal  points  whenever  complete  restocking  occurs.  We  are  concerned 
with  the  probability  distribution  of  the  time  to  absorption,  which  occurs 
whenever  X^(t)  >  s^  for  some  i.  In  principle,  the  steady  state 
probabilities  and  the  expected  tiae  to  absorption  can  be  determined 
analytically  for  these  Markov  processes.  However,  in  aany  practical 
probleas,  the  nuaber  of  possible  states  Bakes  the  transition  matrix  for 
this  problem  prohibitively  large.  For  example,  if  there  are  twenty 
different  Jobs  that  aay  be  encountered  and  there  Is  enough  inventory  so 
that  each  job  aay  be  served  up  to  four  times,  the  total  nuaber  of  different 
states  for  the  vector  N(t)  is  205  -  3.2  Million.  Thus  only  in  very 
small  probleas  will  the  exact  solution  of  the  Markov  process  be  practical. 

Monte  Carlo  Approximation 

An  alternative  approach  for  calculating  the  probability  P{X(t)  &  s)  is 
Monte  Carlo.  This  is  considerably  more  effective  from  a  computational 
standpoint  than  exact  calculation  and  the  nuaber  of  samples  to  obtain  any 
desired  level  of  accuracy  can  be  determined.  For  a  given  job  matrix  J 
and  stock  level  s,  let  us  define  the  indicator  random  variable 

1  if  X(t)  S  s 

Zk" 

0  otherwise , 

corresponding  to  the  kc^  Monte  Carlo  sample.  Then  clearly 

n 

P(X(t)  S  s)  -  Ua  Sn  -  (1/n )l  Zk.  (2.5) 

k-1 

Also,  since  the  variance  of  the  probability  estimate  Sn  is  maximized  when 
the  probability  is  1/2,  we  see  that 

Var(Sn)  -  (l/n)P(X(t)Ss)tl  -  PlX(t)Ss)]  s  l/4n. 

Thus  to  obtain  an  accuracy  of  roughly  two  standard  deviations  (95%)  on  the 
second  decimal  place  of  F(X(t)  s  s),  one  would  need  to  take  n  such  that 

2/<V5n>  at  0.01  or  n  k  10,000. 
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This  number  of  siaulstion  saaples  is  feasible  for  fairly  large  problems, 
and  often  less  accuracy  is  acceptable. 

The  main  difficulty  with  Monte  Carlo  is  that  it  provides  only  a  point 
estimate  and  gives  no  indication  of  how  the  performance  varies  as  a 
function  of  changes  in  stock  level.  This  sensitivity  estimation  is  crucial 
for  inventory  policy  selection  and  optimization.  In  Section  3,  we  will 
derive  a  lower  bound  objective  function  that  is  simpler  to  calculate.  It 
also  provides  the  derivative  information  necessary  for  optimization. 

The  Effects  of  Item  Demand  Dependencies 

In  multi- item  inventory  control  systems,  a  common  assumption  is  that  part 
demands  are  independent  Poisson  arrivals,  with  their  arrival  rates 
determined  by  the  total  item  demands  of  each  type  observed  over  a  period  of 
time.  It  is  interesting  to  consider  the  effects  of  simply  ignoring  the 
interdependencies  between  the  demands  for  items  and  analyzing  the  inventory 
system  as  though  it  had  independent  item  demands.  To  illustrate  this,  let 
us  consider  the  example  in  Figure  2.1.  Let  us  assume  that  jobs  are  Poisson 
arrivals  with  rate  A  -  1,  and  all  jobs  are  equally  likely.  In  this  case, 
the  demand  rates  A^  below  would  be  observed  for  each  of  the  items  i, 
based  on  (2.1). 

Pj  -  1/6.  j  -  1.  2,  3,  4,  5.  6 

Item  1123456789  10 

A£  1.17  0.5  0.5  0.83  0.67  1.67  0.67  0.67  0.83  0.67 

Ve  will  determine  the  probability  of  no  stockout  up  to  time  t  using  the 
independent  item  demand  approach  and  the  true  demand  process  and  compare 
the  results.  In  this  case,  the  problem  is  small  enough  that  an  accurate 
estimate  of  P(X(t)  s  s)  can  easily  be  determined  by  Monte  Carlo. 

The  results  for  the  stock  levels 

s-83354104454 

are  shown  graphically  in  Figure  2.2.  [This  stock  level  was  chosen  to  be 
twice  the  expected  demand  rate  per  unit  {time  for  each  item  type.]  The 
"Exact"  probability  is  P(X(t)  &  s)  evaluated  by  Monte  Carlo.  "Ind  Farts" 
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Exact  Probability  &  Approximations 

Poisson  Demands  Lambda  =  1 

—  Exact  - Ind  Parts  - Parts/ Job 


;novpois  °N  A;iiiqeqoJd 


Elapsed  Time 


1*  the  probability  determined  by  the  product  of  simple  Poisson 
probabilities  that  the  demand  for  part  i  is  less  than  or  equal  to  e^ , 

with  demand  rates  determined  from  (2.1).  The  third  line  in  Figure 

2.2,  labeled  "Parts/Job"  is  based  on  a  bound,  which  we  describe  in  the  next 
section. 

It  is  interesting  to  note  in  Figure  2.2  that  the  probability  of  no 
stockout  calculated  by  the  simple  Poisson  assumption  is  higher  than  the 
Exact  probability  of  no  stockout  for  the  smaller  values  of  t,  while  for 
larger  t  values  it  becomes  less  than  the  Exact  value.  The  errors  in  both 
these  regions  are  clearly  significant.  Thus  the  Interdependencies  in  the 
part  demands  clearly  cannot  be  Ignored  in  many  problems. 

The  "Parts/Job"  bound  has  two  advantages  over  the  "Ind  Parts" 

approximation.  First,  in  the  critical  region  above  90%,  it  is  much  more 

accurate.  Second,  it  always  lies  below  the  true  probability.  Thus,  if 
sufficient  Inventory  is  provided  to  meet  the  service  level  constraint 
calculated  form  the  lower  bound,  we  can  be  certain  that  the  true  service 
level  has  been  met  as  well. 

3.  Obtaining  Bounds 

In  this  section  we  derive  the  "Parts/Job"  bound  in  Figure  2.2.  Since  the 
exact  calculations  of  the  probability  distribution  of  the  time  to  stockout 
are  complex  and  the  sensitivity  to  stock  level  changes  is  even  more 
difficult  to  obtain,  we  seek  an  approximation  for  the  objective  function 
that  makes  it  easy  to  evaluate  and  also  provides  derivative  information. 

Let  us  consider  the  arrival  process  of  repair  jobs  or  demands  for  service, 
which  is  defined  by  the  m  independent  Markov  processes  Nj  (t) .  Ve 
use  the  theory  of  associated  random  variables  to  obtain  a  lower  bound  for 
the  probability  of  completing  all  jobs  up  to  time  t. 

Theorem  1.  If  the  Hj (t)  are  Independent  Markov  processes,  then  we  have 
the  inequalities 

n  n 

1  -^BJl-P(Xl(t)  s  st)l  fc  P(X(t)  5  s)  k^PlX^t)  s  s£). 


-8- 


(3.1) 


Proof:  We  show  that  if  tha  (Nj(t)}  ara  Independent  processes ,  then  for 
each  t,  the  random  variables  (X^ (t) }  are  associated.  Since  the  N,  (t) 
are  independent,  they  are  associated  randoa  variables.  Since  2  0  for 
all  i,j  it  follows  that  the  (X^(t))  are  nondeereasing  functions  of  the 
variables  (Nj(t>).  Therefore,  the  (X^(t))  are  associated  randoa 
variables  themselves.  For  associated  randoa  variables,  the  inequalities  in 
(3.1)  are  known  to  hold.  QED. 

[See  Barlow  and  Proschan  (1975),  pp.  29-34  for  the  developaent  of  the 
results  used  above . ] 

Ve  can  now  focus  our  attention  on  the  calculation  of  the  probability 
distributions  of  the  separate  processes  (X^(t)).  This  in  effect 
substitutes  the  calculation  of  a  one  dimensional  probability  evaluations 
for  one  a  diaensional  probability  evaluation,  which  is  in  general  a  great 
simplification. 


For  one  particular  case  of  interest,  the  evaluation  of  the  individual 
(X^(t))  results  in  a  siaple  closed  fora  solution.  If  the  job  matrix  J 
contains  only  0's  and  l's  ,  that  is,  no  aore  than  one  of  each  part  type 
is  required  per  job  and  the  job  arrivals  are  independent  Poisson  events,  we 
can  express  P(X^(t)  s  s^)  in  closed  fora.  This  is  based  on  the 
observation  that  the  part  demands  theaselves  ara  Poisson  arrivals  in  this 
case.  That  is,  part  i  is  demanded  at  time  t  if  and  only  if  a  job 
arrives  at  time  t,  (a  Poisson  arrival)  and  the  job  is  one  for  which  - 

1.  This  creates  a  Poisson  process  with  rate 

Ai  “  X J  JijPj  • 


Thus  we  have  proved  the  following  result. 

Corollary  1.  If  the  Nj(t)  are  Independent  Poisson  processes  and  the 
aaxiaua  entry  in  J  is  1,  then  we  have 


P(X£(t)  S  al) 


-  jVv  (Att^Al- 

k-0 


(3.2) 


Furthermore ,  by  taking  the  product  of  the  FiX^t)  &  s^ }  and  Integrating 


over  C,  we  obtain  a  lower  bound  for  the  expected  tine  to  first  stockout. 
*1  *n/  k  \  n  . 

'  *  X  •  •  •  l(k1-..kn](l/X)n  (q^i  (3.3) 

kl"°  kn“°  i-1 

where  q^  -  Aj/A  and  k  -  Xjkj. 

Under  the  assumptions  of  Corollary  1,  the  calculation  of  P(X^(t)  s  s^) 
reduces  to  the  simple  "Parts/Job"  method  discussed  previously.  Thus,  in 
this  case,  the  simple  part  fill  approach  provides  a  lower  bound  on  the 
system  performance . 

Calculating  the  Distribution  of  X^(t) 

In  the  general  case,  each  row  in  the  Job  matrix  would  determine  the 
distribution  of  the  number  of  parts  of  type  1  required  when  a  demand 
arrives.  That  is,  we  define  the  random  variable 

Dj(k)  -  the  number  of  items  of  type  I  required  on  the  kth  demand. 

Its  probability  distribution  is  clearly  given  by 

P(D1(k)  -  JjjJ  -  pj,  J  -  1 . m  and  all  k.  (3.4) 

In  the  Poisson  arrival  case,  this  would  mean  that  the  part  demands  are 
Independent  arrivals  with  a  compound  Poisson  distribution.  The  random 
variable  X^(t)  can  be  expressed  as 

M(t) 

Xi(t)  -XDi(k),  (3.5) 

k-1 


where  M(t)  -  Y  Nj(t) . 


In  the  general  case,  we  can  calculate  P(X£(t)  &  s^ )  as  follows.  First 

define 


4 
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(3.6) 


K 

d^K)  -  l  Di(k). 
k-0 


Then  ve  use  the  expansion  formula 


PlX^t)  S  Si)  -  X  PldjOC)  S  Sl)P{M(t)  -  K). 
K-0 


(3.7) 


Then  define 


Fj/KIs^  -  P(di(K)  S  Si),  (3.8) 

that  is,  the  probability  that  part  i  does  not  stock  out  after  K  jobs. 
In  the  general  case  Fj(K|Sj)  is  a  multinomial  probability,  which  must  be 
summed  over  a  particular  set  of  vectors  of  arriving  jobs.  That  is,  ve 
define  the  set  of  vectors 


S1(K|s1)  -  {(nj . i^)  |  X  nj  "  K  *nd  X  njJij 


S  Sj) 


(3.9) 


Then  we  have 


f  K 

F^KISj)  -  X  V»i  ... 
n<Sl(K|si) 


n  (P«)nj  • 

J  J 


(3.10) 


Once  the  F^KIs^)  are  obtained,  P{Xji(t)  S  s^)  is  calculated  from  the 
expansion  formula.  It  should  be  noted  that  although  the  sum  is  infinite, 
the  terms  generally  become  small  quite  quickly.  For  example,  if  the  jobs 
are  Poisson  arrivals,  the  size  of  the  error  from  dropping  the  remaining  terms 
can  be  bounded  by  the  tail  area  of  the  remaining  Poisson  terms  times  the 
last  F^KIsp,  because  F^(K|s^)  is  clearly  decreasing  in  K.  It  is 
also  interesting  to  note  that  truncating  the  expansion  formula  results  in 
obtaining  a  job  completion  probability  that  is  too  small,  so  that  the  error 
from  the  approximation  always  results  in  a  conservative  bound. 


In  the  case  for  which  Jjj  consists  of  only  0’s  and  l's,  F1(K|s1)  can 
be  expressed  more  simply.  As  before,  we  let  q^  equal  the  probability 
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that  part  1  is  demanded  on  a  given  job, 


qi  “  ^  Jij  Pj  "  Ai/A> 

Then 

Fi(K|Si>  “J^i”  (1  *  <Ji)K*“. 


(3.11) 


and  P{Xi(t)  S  s^}  is  given  by  (3.2). 

In  the  general  case,  F^(K|s£)  can  be  calculated  in  principle,  but  it  is  a 
difficult  calculation  for  large  numbers  of  jobs  and  parts.  [It  should  be 
noted,  however,  that  it  is  still  much  simpler  than  calculating  the  joint 
probability  P(X(t)  is).] 

For  the  more  difficult  cases  Monte  Carlo  provides  a  practical  method  for 
determining  F^tKIs^).  In  particular,  one  can  generate  random  sequences  of 
K  jobs  and  simply  observe  the  relative  fraction  of  the  sequences  that  do 
not  produce  stockout.  Since,  this  is  an  estimate  of  a  probability,  the 
variance  of  the  estimate  is  less  than  or  equal  to  l/4n,  where  n  is  the 
number  of  sequences  tested  as  noted  in  (2.5).  This  allows  n  to  be 
selected  to  achieve  the  desired  degree  of  accuracy.  [Note:  In  Section  2, 
Monte  Carlo  was  used  to  obtain  P(X(t)  s  s)  directly,  which  would  be 
quicker  than  evaluating  all  the  F^(K| s^ .  However,  using  Monte  Carlo  to 
evaluate  F^(K| s^)  as  discussed  above  provides  derivative  information 
regarding  the  effects  of  changing  the  stock  level  s^ . ] 

In  Figure  2.2,  the  probability  of  no  stockout  has  been  plotted  as  a 
function  of  elapsed  time  for  the  example  job  matrix  shown  in  Figure  2.1. 
The  "Exact"  and  "Ind  Parts"  lines  were  discussed  at  the  end  of  Section  2. 
The  "Parts/Job"  line  was  obtained  by  evaluating  F^(K| s^)  for  the  given 
job  matrix  and  then  combining  it  with  Poisson  probabilities  P(M(t)-K) 
using  the  expansion  formula.  For  the  time  interval  shown  here  K-6  was 
sufficient  to  obtain  three  decimal  place  accuracy.  Note  that  the 
"Parts/Job"  line  provides  a  pessimistic  (lover)  bound  for  the  exact 
probability  of  no  stockout,  as  Implied  by  Theorem  1.  In  this  case  the 
bound  is  fairly  accurate.  In  addition  since  its  slope  appears  to 
approximate  the  exact  objective  function  fairly  well,  it  should  serve  well 
as  a  surrogate  objective  function.  The  "Parts/Job”  line  in  Figure  2.2  is 
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•Iso  clearly  •  major  Improvement  over  the  "Ind  Parts"  line  In  certain 
regions.  For  example,  when  the  "Ind  Parts"  line  predicts  a  stockout 
probability  of  1%,  the  true  probability  Is  roughly  5%,  while  the 
"Parts/Job"  line  predicts  roughly  6.5*. 

4.  Optimizing  s 

It  Is  important  to  consider  how  F^CKjs^)  varies  with  s^  in  order  to 
solve  optimization  problems.  In  the  case  for  which  contains  only 

0's  and  l's,  it  is  clear,  when  s^  >  Kq^,  that  further  increases  in  s^ 
provide  decreasing  marginal  increases  in  Fi(K|S|).  This  is  because  the 
binomial  distribution  is  unlmodal  with  its  maximum  naar  Kq^ .  However,  in 
the  general  case,  this  monotonicity  does  not  hold.  This  can  be  illustrated 


intuitively  by  considering  the  row 
probabilities  shown  below. 

Jil . Jim 

and  associated 

#  part  1  required  1 

2 

3 

4 

Pj  0.01 

0.01 

0.01 

0.97 

Clearly  as  s^  passes  through  each  multiple  of  four,  there  will  be  a 
substantial  jump  in  the  probability  of  no  stockout.  Thus  the  property  of 
decreasing  marginal  improvements  holds  only  in  the  simpler  case. 
Furthermore,  this  same  type  of  argument  can  be  used  to  show  that  the 
original  objective  function  P{X(t)  s  s)  is  not  concave  in  s  either. 

In  choosing  the  stock  level  s,  one  would  typically  consider  an 
optimization  problem  of  the  following  form 

min  T  hjs^,  subject  to  P(X(t)  s  s)  fc  a.  (4.1) 

l.i>  1 

This  problem  has  Lagrangian 

*<•)  -  I  hlBi  -  pP<X(t)  s  a).  (4.2) 


This  Lagrangian  is  clearly  analogous  to  the  one  that  results  from  the 
problem 
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■ax  P(X(t)  s  a)  subject  to  Y  hjSj  5  b. 
“i>  T 


(4.3) 


These  problems  are  typically  not  computationally  feasible,  because 
of  the  difficulty  in  evaluating  the  objective  function.  Furthermore,  as 
noted  previously,  the  objective  functions  may  not  be  concave  in  a. 


An  alternative  optimization  approach  for  (4.3),  or  with  the  appropriate 
modification  (4.1),  is  to  consider  meeting  the  lover  bound  on  service  level 
from  Theorem  1  at  minimum  total  inventory  cost.  That  is,  we  consider  the 
problem 


■ax  £  log  P(Xi(t)  S  si)  subject  to  £  h^s^  s  b, 


(4.4) 


which  has  Lagrangian 


iKs)  -  X  log  PlXi(t)  s  s, }  -  /iX  hisi 
i  1  i 


(4.5) 


Problem  (4.5)  has  the  major  advantage  of  being  additively  separable.  In 
the  simpler  case,  for  which  Jjj  -0,1,  the  Individual  functions 
log(P{X1(t)  S  8ll)  -  ph^s^  are  concave  in  s^.  Let  us  see  why  this  is 

true.  First,  as  noted  above,  all  the  functions  Fj(K|Sj)  experience 
decreasing  marginal  improvements  as  s^  is  increased  beyond  s^  >  KAj/>. 
Since  PtX^t)  5  s^)  is  a  weighted  sum  over  K  of  the  F1(K|si),  where 
the  weights  are  invariant  with  s^ ,  P(X^(t)  s  s^}  will  experience 

decreasing  marginal  improvements  as  well.  Thus  a  marginal  analysis 
approach  of  allocating  additional  stock  s^  to  the  part  that  has  the  best 
improvement  per  unit  cost  will  lead  to  a  unique  optimum  when  Jjj  S  1  for 
all  i.J. 

Example  Optimization 

Let  consider  the  solution  of  a  problem  of  the  form  (4.5)  to  illustrate  the 
solution  technique  and  the  accuracy  of  the  approximation.  For  this 
example,  we  will  not  use  the  job  matrix  in  Figure  2.1,  because  it  has 
entries  greater  than  one,  which  would  cause  the  individual  objective 
functions  in  (4.5)  to  have  multiple  local  optima.  Instead,  we  will 
consider  the  job  matrix  shown  below  with  all  Jjj  s  1. 
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Jobs 

111111 
111110 
11110  0 
parts  1110  0  0 

1  1  0  0  0  0 
1  0  0  0  0  0 

The  holding  costs  {h^}  per  unit  will  be  1,1,2, 2,3, 3,  respectively.  Jobs 
are  as suae d  to  be  Poisson  arrivals  with  rate  X  -  1  and  each  job  type  is 
equally  likely  to  occur.  Thus  the  individual  part  demand  rates  per 

unit  time  are,  respectively 

1,  5/6,  2/3,  1/2,  1/3,  1/6. 

We  will  take  t  -  6  and  solve  problem  (4.5),  with  P(X^(6)  s  s^}  given  by 
(3.2).  Instead  of  choosing  a  particular  budget  constraint,  we  will  develop 
the  exchange  curve  for  probability  of  no  stockout  before  t  -  6  for  all 
inventory  budgets  greater  than  50. 

First  we  need  to  determine  the  form  of  the  optimal  stocking  policy  for 
(4.5)  for  this  problem.  Since  we  are  developing  an  exchange  curve,  what  we 
really  need  is  the  optimal  order  in  which  to  stock  the  parts  as  the  budget 
is  Increased.  This  can  be  determined  as  follows.  Consider  the  value  of 
p,  call  it  p(s^),  such  that  we  become  indifferent  between  the  stock 
level  8^  and  the  stock  level  Sj+1  for  part  i.  If  we  consider  A^(s^) 

-  0  for  the  Lagranglan  in  (4.5),  we  obtain 

p(sA)  -  [log  P(X1(6)  S  sj+l)  -  log  PlXi(6)  s  si)]/hi  (4.6) 

for  each  part  1  and  stock  level  s^ .  Because  of  the  nature  of  the  Poisson 
probability  in  (3.2),  each  p(st)  is  strictly  decreasing  in  «£.  We  can 
use  the  right  hand  side  of  (4.6)  to  rank  the  stock  levels  in  terms  of  the 
value  of  p  chat  will  cause  them  to  first  be  activated.  Clearly  as  p 
decreases  we  pass  through  the  various  stock  levels  of  all  parts  in  the 
optimal  order,  since  the  same  p  serves  for  all  parts  simultaneously. 

This  ordering  of  right-hand  sides  of  (4.6)  across  all  part  types  and  stock 
levels  then  gives  the  optimal  order  for  successively  increasing  the  stock 
levels  in  solving  (4.5). 
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The  resulting  exchange  curve  for  this  problea  is  s* o'"  in  Figure  4.1.  This 
graph  vas  constructed  from  a  few  saaple  stock  levels,  as  shown  below, 
although  all  stock  levels  were  deterained  in  sequence  for  the  problem. 


Inventory 

Cost 

{Vi 

*1 

*2 

Stock 

Levels 

*3  s4  1 

b5 

*6 

50 

9 

8 

5 

4 

3 

2 

56 

10 

9 

6 

5 

3 

2 

63 

10 

9 

7 

6 

4 

2 

68 

11 

10 

7 

6 

4 

3 

79 

12 

11 

9 

7 

5 

3 

88 

13 

11 

9 

8 

5 

4 

98 

14 

13 

10 

9 

7 

4 

105 

15 

14 

11 

9 

7 

5 

The  exchange  curve  in  Figure  4.1  shows  the  probability  of  no  stockout  that 
was  achieved  by  each  of  these  stock  levels  for  the  objective  function  in 
(4-3). 

To  test  the  accuracy  of  the  bound  in  Theorem  1,  the  actual  job  completion 
probability  P(X(6)  s  s)  was  also  deterained  for  each  of  the  stock  levels 
s  shown  above  by  Monte  Carlo.  This  is  also  plotted  on  Figure  4.1.  [N.B. 

This  actual  probability  calculation  used  the  stock  levels  s  shown  above, 
which  are  not  the  optima  for  the  actual  objective  function! ]  In  this 
example  the  difference  between  the  bound  and  the  actual  objective  function 
is  considerably  larger  than  it  was  for  the  job  matrix  in  Figure  2.1.  This 
is  probably  because  the  job  matrix  in  this  section  has  much  more 
commonality  or  linking  of  parts  across  jobs,  making  the  independent 
treatment  of  parts  a  poorer  assumption. 

However,  from  an  optimization  standpoint,  we  may  read  the  probability  of 
no  stockout  from  the  actual  curve  to  determine  the  performance  that  is 
achieved  for  any  given  inventory  investment.  The  error  introduced  by  our 
approximation  of  using  the  bound  instead  of  the  true  objective  function  is 
really  the  difference  between  the  performance  that  the  Inventory  investment 
in  Figure  4.1  achieves  and  the  performance  that  could  have  been  achieved  by 
optimizing  the  stock  level  for  the  actual  objective  function  subject  to  the 


£4 
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Ill 


Cost 


•a»  budget.  Because  of  the  difficulties  in  optimizing  the  actual 
objective  function,  there  is  no  way  to  determine  precisely  how  large  this 
error  is.  However,  for  any  budget  level,  a  corresponding  stock  level  and 
probability  of  no  stockout  obtained  fron  Figure  4.1,  we  can  be  certain  of 
achieving  the  performance  level  shown  with  that  inventory  investment. 

5.  Conclusion 

This  paper  considers  the  problem  of  selecting  stock  levels  for  multi- items 
to  meet  a  combined  objective  of  job  completion  rate  or  probability  of 
filling  all  orders  when  demands  for  parts  are  interrelated.  As  noted  by 
previous  authors,  this  is  both  a  difficult  and  important  problem  in  many 
Inventory  systems.  Ve  believe  that  exact  solutions  will  be  computationally 
unfeasible  for  the  optimization  problems  resulting  from  many  applications. 
We  have  developed  a  lower  bound  function  for  the.  service  level  objective 
which  is  simpler  to  calculate.  In  the  case  of  job  matrices  with  all  0's 
and  l's  this  lover  bound  objective  function  is  also  relatively  easy  to 
optimize.  The  stock  levels  resulting  from  optimization  of  the  lover  bound 
can  then  be  evaluated  for  the  true  objective  function  using  Monte  Carlo. 
This  allows  an  exchange  curve  to  be  derived  for  the  true  objective 
function,  using  reasonable  (although  generally  not  optimal)  stock  levels. 
This  combination  of  optimizing  a  surrogate  objective  function  and 
evaluating  the  true  objective  function  by  Monte  Carlo  appears  to  be  an 
attractive  approach  for  a  class  of  problems  which  would  not  be  solvable 
otherwise . 

Clearly  there  are  important  classes  of  problems  that  are  beyond  the  scope 
of  the  solution  methods  described  in  this  paper.  Therefore,  ve  believe 
that  this  is  fruitful  area  for  further  research,  both  In  the  context  of  Job 
completion  rates  for  repair  service  systems  and  for  order  fill  rates  in 
manufacturing  systems  with  parts  commonality. 
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